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We discuss a simple, experimentally feasible scheme, which 
elucidates the principles of controlling ("engineering") the 
reservoir spectrum and the spectral broadening incurred by 
repeated measurements. This control can yield either the in- 
hibition (Zeno effect) or the acceleration (anti-Zeno effect) of 
the quasi-exponential decay of the observed state by means 
of frequent measurements. In the discussed scheme, a pho- 
ton is bouncing back and forth between two perfect mirrors, 
each time passing a polarization rotator. The horizontal and 
vertical polarizations can be viewed as analogs of an excited 
and a ground state of a two level system (TLS). A polariza- 
tion beam splitter and an absorber for the vertically polarized 
photon are inserted between the mirrors, and effect measure- 
ments of the polarization. The polarization angle acquired in 
the electrooptic polarization rotator can fluctuate randomly, 
e.g., via noisy modulation. In the absence of an absorber the 
polarization randomization corresponds to TLS decay into an 
infinite-temperature reservoir. The non-Markovian nature of 
the decay stems from the many round-trips required for the 
randomization. We consider the influence of the polarization 
measurements by the absorber on this non-Markovian decay, 
and develop a theory of the Zeno and anti-Zeno effects in this 
system. 

PACS numbers: 03.65. Bz, 42.50.-p, 03.67.-a, 42. 25. Hz 



I. INTRODUCTION 

The quantum Zeno effect (QZE) is the striking predic- 
tion that the decay of any unstable quantum state can be 
inhibited by sufficiently frequent observations (measure- 
ments) The QZE has been experimentally tested 
H and primarily analyzed for two coupled states (^-0] 
(with few exceptions |^,^|). Yet the consensus opinion 
has upheld the QZE as a general feature of quantum 
mechanics which should lead, e. g., to the inhibition of 
radioactive or radiative decay [jTo| JlTf] . The claim of QZE 
generality has rested on the assumption that successive 
observations can, in principle, always be made at time 
intervals too short for the system to change appreciably. 
However, this assumption and the generality of the QZE 
have scarcely been investigated. 

We have now shown fl2|j that this assumption is ba- 
sically incorrect and that the QZE does not hold gen- 
erally, but only in a restricted class of systems. The 
main implications of our theory are: (i) The QZE is 
principally unattainable in radiative or radioactive de- 
cay, because the required measurement rates may cause 



the system to disintegrate (via the production of new 
particles), (ii) Decay acceleration by frequent measure- 
ments (the anti-Zeno effect - AZE) |l2|-|l5fl is possible 
for essentially any decay process, and is thus much more 
ubiquitous than its inhibition (the QZE). These findings 
stem from the universal result |T^ ] , whereby the modifi- 
cation of the decay rate by frequent measurements is de- 
termined by the convolution (overlap) of two functions: 
(a) the measurement-induced spectral broadening (en- 
ergy spread), which is proportional to the rate of mea- 
surements, in accordance with the time-energy uncer- 
tainty relation; (b) the spectrum of the reservoir (bath) 
to which the decaying state is coupled. The QZE or 
AZE correspond to the measurement-induced spread be- 
ing much broader or narrower than the reservoir spectral 
width, respectively. The non-Markovian nature of any 
physical decay process, associated with the finite spectral 
width (and, correspondingly, nonzero memory time) of 
the reservoir, is the essential property allowing the mod- 
ification of the decay by means of frequent measurements, 
be it the QZE or the AZE. 

The universal formula jl2| was obtained on the basis 
of the projection postulate. In reality there can be many 
different measurement schemes, which can be classified as 
direct, when the initial state itself is measured |JJ|,@,§|, 
and indirect |11 T^,|l^l, when the final state(s) are mea- 
sured. Whereas direct measurements should be nonde- 
structive, i.e., conserving the population of the measured 
state, indirect measurements can be either nondestruc- 
tive or destructive. This difference can affect the dynam- 
ics for large times, when the initial-state population sig- 
nificantly differs from 1. However, for short times, when 
this population is close to one, formula obtained in [T^ ] 
should hold for the both types of measurements. 

In this paper we discuss a simple, experimentally feasi- 
ble scheme, which elucidates the principles of controlling 
( "engineering" ) the reservoir spectrum and the spectral 
broadening incurred by repeated measurements. This 
control can yield either the inhibition (QZE) or the accel- 
eration (AZE) of the quasi-exponential (non-Markovian) 
decay of the observed state by means of frequent mea- 
surements. In the present scheme, the pertinent observ- 
able is photon polarization, constituting the optical ana- 
log of a two- level system [jl6). The photon is injected into 
the setup via a fast gate, which first reflects and then 
transmits horizontal polarization. The injected photon 
is bouncing back and forth between two perfect mirrors, 
each time passing a polarization rotator (Fig. [J18| . A 
polarization beam splitter (PBS), which reflects a verti- 
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cally polarized photon and transmits a horizontally po- 
larized one, as well as an absorber, are inserted between 
the mirrors. If the path of the vertically polarized photon 
is completely blocked by the absorber, then it realizes a 
discrete measurement (projection) at each passage: the 
photon is cither lost or its state is projected onto the 
horizontal polarization. This situation has been used to 
demonstrate an interaction-free measurement |l9| --plf. If, 
on the other hand, the absorber is partially transparent 
to vertically polarized photons, then it realizes an imper- 
fect measurement [ p2|p5| ]. 

The polarization angle acquired in the electrooptic po- 
larization rotator can fluctuate randomly, e.g., via noisy 
modulation of a Pockcls cell. In the absence of the ab- 
sorber, the polarization, after many round-trips, then be- 
comes random and the probability of finding any par- 
ticular polarization tends to 1/2. Taking the horizon- 
tal and vertical polarizations as analogs of an excited 
and a ground state of a two level system, this polariza- 
tion randomization corresponds to decay into an infinite- 
temperature reservoir. The non-Markovian nature of the 
decay stems from the many round-trips required for the 
polarization to change randomly, which implies a long 
memory time. Our goal is to consider the influence of 
the polarization measurements by the absorber on this 
non-Markovian decay, and derive the conditions of the 
QZE and AZE in this system. It should be noted that 
if the source of injected photons is a laser, governed by 
quasiclassical photon statistics, then the results of polar- 
ization decay can be interpreted classically. However, the 
concept of the Zeno effect is valid in classical electromag- 
netism, as pointed out by Peres [fTo] . 

In Sec. |l| we present the physical model and its fea- 
tures in the limits of perfect absorption or fixed (rather 
than random) rotation angles. In Sec. Ill the master 



equations for its general analysis are derived (any ab- 
sorption and rotation). Sec. |V| is devoted to continuous 
dephasing, i.e., the limit of small, highly correlated ran- 
dom rotations and weak absorption (ineffective measure- 
ments) . Sections [v| and VI deal with discrete dephasing 
(namely, correlated and anticorrelated phase jumps) and 
arbitrary absorption (effective and ineffective measure- 
ments) . Conditions for the QZE and AZE are derived. In 
particular, in Sec. [y|a general theory is developed for the 
most interesting case of small rotation angles, whereas in 
Sec. VI a simple non-Markovian model for random ro- 



tations of an arbitrary size is studied analyti cally and 
numerically. The conclusions are given in Sec. VII . Ap- 
pendices [A], |b] and ^| give the details of the general analy- 
sis, the small and the arbitrary-size phase-jump analysis, 
respectively. 



II. MODEL DESCRIPTION 



A. The setup 

Consider the setup in Fig. [I] pl|. A horizontally po- 
larized photon, denoted as \h), enters the setup via a fast 
gate (not shown) , which changes from being totally reflec- 
tive to totally transparent to \h) on a ns scale. The polar- 
ization rotator between the two highly reflecting mirrors 
causes fast (ns-scale) rotation of the photon polarization 
by means of a Pockels cell or another electrooptic ele- 
ment. The \h) photon, transmitted by the PBS, bounces 
between the mirrors, while a vertically polarized photon, 
denoted as \v), which is reflected by the PBS, is blocked 
by an absorber, which can be made partially transparent 
with transmissivity 9. (Alternatively, one can use a per- 
fect absorber and a PBS which is partially transparent 
for the vertically polarized photon \v).) 
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FIG. 1. Setup for controlling the polarization decay of a 
photon bouncing between the mirrors. Measurements are ef- 
fected by a polarization beam-splitter (PBS) and an absorber 
with variable transparency 9. The "reservoir" into which the 
polarization decays is realized upon modulating a Pockels cell 
(which rotates the polarization by Aip) by a field with con- 
trollable noise properties. 



B. QZE for constant rotation angle and complete 
absorption 

When the path of \v) is open, the probability of finding 
the photon with horizontal polarization after n round 
trips is 



Ph{n) = cos 2 (nAip) , 



(1) 



assuming a fixed Aip, the rotation angle of the polarizer 
at each passage: the polarization oscillates between the 
vertical and horizontal states, analogously to a Rabi oscil- 
lation. If, on the other hand, the path of \v) is completely 
blocked by a perfect absorber, this absorber realizes an 
impulsive measurement (IM), i.e. a projection, at each 
passage: the photon is either lost or its state is projected 
onto the \h) state. With the completely blocked \v) path 
the probability of finding \h) decays as 



Ph(n) 



cos 



2ii 



(2) 
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where n is the number of the round-trips. With Aip> -C 
1 the effective exponential decay is much slower than the 
Rabi oscillation. This constitutes an example of the QZE 
by IMs |4U^ ]l^ , |l3| , ^lf , which has been introduced as a 
demonstration of interaction- free measurements |T^-^J . 



C. Arbitrary absorption and rotation 

To extend the above results to the case of partial ab- 
sorption and arbitrary rotation angle Aipk (fc = 1, . . .), 
we write down the obvious dynamical equations for the 
horizontal and vertical components of the photon field 
envelope (wave function), eh and e v , respectively, 



e h = -fj,(t)e v , 

r 0) 

e v = fi(t)e h — e„. 



(3) 



Here /j,(t) = ip(t), cp(t) is the total angle accumulated 
by consecutive Atpk between the polarization vector and 
the vertical axis of the polarizer, and T(t) is the photon 
absorption rate in each passage. The rates of polarization 
rotation and absorption, fj,(t) and T(t), are time-periodic 
functions, whose period is the round-trip time r r . They 
do not temporally overlap in the setup of Fig. |l|, i.e., 
there exists t\ (t\ < r r ) such that n(t) vanishes for t\ < 
t < r r while r(t) vanishes for < t < T\. Equation (|^) is 
equivalent to the Schrodinger equation for an open two- 
level system (TLS) interacting with afield (J,(t). 
The general solution of Eq. (0) at t = nr r is 



where 



n 

,fe=i 
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9 = cxp 



T(t)dt 



(4) 



(5) 



is the amplitude transmission coefficient of the absorber. 
Within the above constraints /i(t) and T(t) are arbitrary. 
Most of the subsequent analysis will be dedicated to ran- 
dom fi(t) caused by noisy modulation of the electro-optic 
rotator. 



D. QZE for fixed rotation angle and incomplete 
absorption 

Assuming again a fixed Aip (all Aipi- being equal), one 
obtains from (Q) the following time-dependent probabil- 
ity for the photon to keep its initial horizontal polariza- 
tion (e h (0) = 1, e„(0) = 0) 



Ph(t = nr r ) = e 2 h (t) 

= {[A?(cosA^ - A 2 ) + A£ (Ai - cos A^)]/D} 2 , 



where D = ^(l + 6») 2 cos 2 Aip - 49 and Ai, 2 = [(1 + 
9)cosA<p±D]/2. For complete transparency, 9=1, Eq. 
(0) reduces to the result ([j]) and for complete absorption, 
9 = 0, to the result (§). 

For small phase jumps Atp and sufficiently strong ab- 
sorption, (A(p) 2 <C (1 — 9) 2 , the photon polarization de- 
cay is approximately exponential 



P h (t = UT r ) = exp 



(l + g)(Ay>) ; 
(1 - 0) Tl 



(7) 



It is intuitively clear that the effective rate of measure- 
ments v increases with the quantity 1 — 9, which plays 
the role of the effectiveness of the measurements. This 
will be rigorously confirmed in Sec. Correspondingly, 
the decay rate in Eq. (]?]) decreases with the increase of 
1 — 9, thus demonstrating the QZE in the case of fixed 
Aip. 



III. GENERAL ANALYSIS OF POLARIZATION 
DEPHASING: MASTER EQUATIONS 

In the ensuing analysis we assume that the rotation 
angle Aip is not fixed but fluctuates randomly, e.g., due 
to the modulation of the Pockels cell rotator by a noisy 
control field. Thus, after many round trips the polariza- 
tion of the photon becomes random and the probability of 
finding any particular polarization tends to 1/2. Taking 
the horizontal and vertical polarizations as analogs of the 
excited and ground states of a two level system (TLS), 
this polarization randomization corresponds to decay into 
an infinite-temperature reservoir. 

In what follows, we shall write down the most general 
equations describing the influence of polarization projec- 
tion measurements on such decay. The polarization ten- 
sor (or the density matrix for the TLS) obeys the equa- 
tion 



Q = [A(t) + Cn(t)]Q. 



(8) 



Here Q = (Ph, P v ,v)^ , P v (t = nr r ) = e 2 (t) is the prob- 
ability for the photon to have the vertical polarization, 
u = 2e v Eh is the coherence, in TLS terms, or the first 
Stokes parameter fp4|, and 



A(t) = 
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(9) 



Equation (||) can be reduced to an equation for the aver- 
age quantity Q(t), involving an expansion in cumulants 
of C/J,(t) p5| , p6| . As shown in Appendix by trun- 
cating the cumulant expansion at the second order, one 
obtains the following non-Markovian differential master 
equations for the average polarization probabilities, 



-R(t)P h + R(t)P v , 



(6) 



dPh 
dt 

dP - 

^ = R(t)P h -[R(t) + T(t)}P v . 



(10) 
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The rate that governs the polarization change in (|icj) is 
given by the integral 

R(t) = 2[ dt'k(t,t')9(t,t'), (11) 
Jo 



whose integrand is the product of 



(12) 



the correlation function of the random rotation rate n(t), 
and of 



0(t,t') 



exp 



1 



r(ti)dti 



(13) 



which is related to the polarizer transparency (measure- 
ment effectiveness, as discussed below). We assume that 
the average of /i(t) vanishes, /j,(t) = (no systematic shift 

To obtain the validity condition of the above master 
equations, one should consider higher-order terms in the 
cumulant expansion. As shown in Appendix |X| for the 
case of continuous noise, the comparison of the second 
and fourth cumulants implies that Eqs. ( [To| ) hold under 
the condition 



(14) 



where Tr is the reciprocal correlation time of fx(t). Below 
we assume that criterion ( pT] ) holds also for the case of 
discrete noise. 



IV. CONTINUOUS DEPHASING 

In this section we consider the case when Ph,v(t) vary 
slowly on the time scale of several round trips, which 
allows one to describe them by continuous functions of 
time on the coarse-grained time scale (a more general 
case will be discussed in Sec. |v|). In the present case the 
phase jumps 



dtn(t) 



(15) 



(n-l)i 



are small, B 2 = ((Aip n ) 2 ) <C 1, and highly correlated, 
Atp n w A(/3„_i, whereas the measurements are highly 
non-effective, 8 « 1. Then T{t) can be substituted by 



r = 



'r Jo 



dtT(t), 



(16) 



whereas n(t) can be considered as a continuous, station- 
ary random process, implying k(t,t') = k(t — t'). As a 
result, now in Eqs. ( |l0| ) 



R(t) = 2 [ dt'k(t') exp 
Jo 



Tot' 
2 



(17) 



For t , where Tr is the characteristic decay rate 

of k(t), Eqs. ( |l0| ) become Markovian (although they ac- 
count for non-Markovian phase fluctuations associated 
with fi(t)) 



dP h 

-f = -RP h + RP v , 
dt 



dP v 
dt 



= RP h -{R + T )P v , 



with constant rate 
R = 2 



J dtk(t)cxp^ — ^- 



(18) 



(19) 



The last factor in the integrand here expresses the decay 
law of the first Stokes parameter, f(t) = u(t)/u(0), due 
to measurement-induced relaxation. This follows from 
Eqs. (|) [with fj,(t) = 0] and ©. 

The solution of Eqs. ([l^) with the initial conditions 



yields 
Ph{t) 



P h (0) = l, il(0) = 0, (20) 



, -(R+T a /2)t f CQshSt + £0 sinhSt | (2] • 



where S = ^R 2 + T 2 /4:. 

In the absence of measurements, Tq = 0, ( ^l|) yields 



P h (t) = |(1 + . 



(22) 



where Rq = 2 dtk(t). By contrast, the measurement- 
affected polarization decays, assuming that Tq ^> R, as 



hit) 



-Hi 



(23) 



so that R is indeed the measurement-modified decay rate. 

The expression for the decay rate (19) can be recast in 
the same form as the universal result in Ref . , 



Here 



R = 2ir 



G(uj)F(uj)duj. 



dtk(t) cos Lot 



(24) 



(25) 



is the random-field intensity spectrum, which can be con- 
sidered as the spectrum of the infinite-temperature reser- 
voir. The other factor in (E4h, 



1 



F(u) = - 



r /2 



^(r /2) 2 + ^ 2: 



(26) 



is the Fourier transform of the measurement-induced de- 
cay law of the first Stokes parameter. The width of F(oj) 
has the meaning of the effective rate of measurements v 
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fl2"| |. The definition of v introduced in |T^] becomes in 
the present case of a degenerate TLS 



(27) 



From Eq. (|2q) one obtains v = To/2. 

Equation (p 4) allows a graphical interpretation of the 
QZE In particular, Eq. (§|) shows that the 

QZE occurs when the reservoir spectrum is peaked around 
u> = 0. However if the spectral peak of the reservoir is 
detuned from the resonance frequency uj a of the TLS (here 
Lu a = 0), one can obtain the quantum anti-Zeno effect 
(AZE), i.e., an increase of the decay rate R wit h the 
effective measurement rate v, as illustrated in Sec. VA. 



V. DISCRETE DEPHASING (PHASE JUMPS) 



One recognizes in this expression the analog of Eq. ( |T9| ) 
for the discrete case, on realizing that 9 n (n > 0) is the 
discrete measurement-induced relaxation function of the 
first Stokes coefficient. 



B. Correlated and anticorrelated discrete jumps: 
QZE and AZE 

The decay rate R ( |3l| ) is independent of the measure- 
ments for uncorrelated jumps, i.e., Markovian phase fluc- 
tuations (7 = 0). However, R as a function of the effec- 
tiveness of measurements 1 — 9 decreases for correlated 
phase jumps (7 > 0), thus demonstrating the QZE, and 
increases for anticorrelated jumps (7 < 0), demonstrat- 
ing the AZE (Fig. 0) . How can we interpret these results? 



A. Small jumps: General analysis 

Here we allow for any degree of correlation between 
consecutive phase jumps, as well as for arbitrary absorp- 
tion per passage. In this subsection a general theory is 
developed for the case of sufficiently small phase jumps. 

We consider the sequence of polarization rotations by 
the angles Aip n to be a discrete-time random process 
with the correlation function K nm — {A(p n Aip m ) . Typ- 
ically, the random process Atp n is stationary, yielding 
K nm — K n ^ m = K rn ^ n . The general analysis is given in 
Appendix |b[ Here we present the simple case 



(28) 



where 7 (— 1 < 7 < 1) is the correlation degree between 
two successive jumps EfJ: Aip n w A(p n+ i for 7 w 1, 
Aip n and A(p n+ i tend to have opposite signs for 7 < 
and are statistically independent for 7 = 0. In thi s case 
[Eq. (p8|)1 the correlation time is given by [cf. ( |f3l|) 1 



TV 



7 



(29) 



We start from the master equations (10), which are 
applicable both for continuous and discrete evolution. As 
shown in Appendix pi under the conditions 



n > 



\j\6 



1 - ( 7 0) 2 



B 2 «(l- 7 )(l-70), 
Eqs. ( |Toj ) yield solution Eq. (pl|), where now 

R=- t — , t = nr t . 

1 — 7# r r 

For the general case 



R = 



1 



R n0 H 



(30a) 
(30b) 

(31) 
(32) 




FIG. 2. Decay rate R dependence on measurement effec- 
tiveness 1 — 8 for different degrees of correlation between con- 
secutive phase jumps, according to Eqs. ( fjij ) and ( |33| ) Here 
B — 0.1, r r = 0.07. For the curves from top to bottom 
7 = 0.7 (correlation leading to QZE), 0, —0.9 (anticorrelation 
leading to AZE) . 

The first step towards gaining insight into the decay 
rate ( |3l| ) is to realize that it can be recast in the form 
(Appendix |b|) 



R = 2tt 



7t/t t 



7r/r r 



dujG{uj)F(uj), 



where 



D 2 



1 



•7 



27rr r 1 + 7 2 — 27 COS UJT T 
is the reservoir spectrum, and 

tv 1-e 2 



F(u) 



2tt 1 + 9 2 - 29 cos lot, 



(33) 



(34) 



(35) 



is the Fourier transform of the measurement-induced re- 
laxation function of the first Stokes parameter [see Eq. 
(p9j)]- The form ( pS3| ) is essentially the same as (|24|) or 
the universal form of measurement-affected decay rates 
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given in Ref. |i~2|. What makes ( |33| ) distinct is that we whereas for highly anticorrelated jumps (7 f» —1), it is 

are now in a finite frequency domain, — 7r/r r < u> < 7r/r r , peaked near lu — ±7r/r r , 

in contrast to Eq. (p4|). This is related to the fact that ^ 

polarization evolution studied here is discrete in time. In G(u) ~ — — - (7r ± lot <C 1) (41) 

what follows, we shall first separately analyze F(w) and ttt 2 r^ 2 + (k/t t ± w) 2 
G(u>), and then graphically deduce the QZE and AZE 

from their convolution (1331). " with r R = ( : + 7)Av- 



_Z. Properties of F(ui) 



3. Graphical analysis of the decay rate 



The Fourier-transformed measurement-induced relax- 
ation function of the first Stokes parameter (|35|) is nor- 
malized to one. It has one maximum, F max (uj = 0) = 
(r r /2vr)(l + ff)/(l - 9) located at lo = 0. It is mini- 
mal at the borders lo = ±7r/r r of the frequency domain, 
-Fmin = (r r /27r)(l — 9)/(l + 9). In the particular case of 
the ideal (fully effective) projective measurements (IMs), 
F(lo) is constant, 



F(lo) 



2tt 



(9 = 0), 



(36) 



whereas for non-effective (unreliable) measurements 
1), F(ui) is a narrow peak [cf. Eq. (26)], 



F{u) 



1 



r /2 



7r(r /2)2 +w 2 



(|w|r r « 1). 



(37) 



Here we took into account that now (1 — 9) ~ Tq^/2 <C 1 
[cf. Eqs. (|) and ©]. 

Inserting Eq. (|35| ) into (^7|), one obtains that the ef- 
fective rate of measurements is 



2(1 



1 + 9 T T 



(38) 



For any allowed value of 9, we have v ~ (1 — 9)/t T) which 
formally confirms the intuitive interpretation of 1 — 9 as 
the effectiveness of measurements. More specifically, for 
IMs, v = 2/t t , which differs by a factor of 2 from the 
real rate of measurements (note that the definition of v 
is meaningful only with an accuracy up to a factor of the 
order of one). For low-effectiveness (highly unreliable) 
measurements (9 ~ 1), we have v = (1 — 9)/t t = Tq/2. 



2. Properties of G(ui) 

The reservoir-coupling spectrum (^) is constant in the 
Markovian case G(u>), 7 = 0, 



G(lo) 



B 2 
2ttt t 



(7 = 0). 



(39) 



By contrast, G(u>) is mainly concentrated near u> = for 
highly correlated jumps (7 « 1), 



G(lo) 



b 2 r R 

7TT-2 T 2 , + L 



(Mr r « 1), (40) 



Since the decay rate Eq. (|33| ) is determined by the 
overlap of F(lu) and G(u>), the above results allow a 
graphical interpretation of the dependence of R on the 
effective measurement rate (Fig. ||). In Fig. ||(a), G(lo) 
is flat (7 = 0) and the convolution ([33]) is proportional 

to the integral of F(uj), JlfjL oIluF(lu) = 1. As a re- 
sult, R is independent of the shape of F(uA (i.e., of the 
effectiveness of measurements). In Fig. 0(b) (7 > 0), 
G(lu) is peaked at u> = and hence the convolution R 
is determined by the portion of F(u) inside the width of 
G(lo). This implies a reduction of R with the broadening 
of F(u>), i.e., the QZE. The opposite is true in Fig. |^(c), 
due to the fact that the reservoir peaks are detuned from 
the TLS frequency Lo a = 0. Now the convolution R is 
sensitive to the wings of F(lo), which rise with the in- 
crease of F(w). As a result, the AZE occurs (increase of 
R with the broadening of F(ui), i.e., with the increase of 
the measurement effectiveness). In Figs. ||(b), ||(c), the 
limit of the flat F(ui) (i.e., IMs) implies that the integral 
([33|) is independent of the shape of G(oj), resulting in the 
same value of R as in the Markovian case, Fig. ||(a). 

The decay rate R corresponding to Eq. (0) can be 
recast as 



R = 2(A ¥ >/r r )>, 



(42) 



with the effective rate of measurements v given by Eq. 
(|38|). The above expression for R has the familiar 
form characteristic of the QZE ]l2| ]. The coefficient 
(Aip/r r ) 2 = p, 2 is the squared average of the coupling 
amplitude over the round trip. 



VI. EXACTLY SOLVABLE MODELS OF 
RANDOM PHASE JUMPS 

We now consider simple models allowing an exact solu- 
tion in the cases of absent measurements (free evolution) 
and perfect measurements (IMs)]. 

In the case of free evolution the probability that the 
photon is found with the horizontal polarization is given 
by 

P h (n) = (cos 2 <p n ) = 1 + iRe(e 2 ^"), (43) 

where ip n = 5Zfe=i A(pk. and the angle brackets denote 
the averaging. 
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FIG. 3. Overlap of the reservoir spectrum G(u>) in ( |j^ ) and 
the repeated- measurement broadening F(w) in (|35|), deter- 



mining the convolution integral of the decay rate 7? in (B5). 
Solid lines: G{uS) with (a) 7 = (uncorrelated jumps), (b) 
7 = 0.7 (correlated jumps - QZE), (c) 7 = —0.7 (anticorre- 
lated jumps - AZE). Dashed lines: F{oj) with 8 — (perfect 
projections). Dotted lines: F(ui) with 8 — 0.9 (ineffective 
measurements). Here B = 0.1, r r = 0.07. 



A. Markovian limit (uncorrelated jumps) 

First, let us assume that the rotation angles Aip during 
the round-trips are independent of each other. We can 
then write Eq. (EJ3) in the form 



x exp 



2i ^2 A( Pk 



fe=i 



dAipi . . . dAip Ti 



(44) 



Assuming the probability distribution of each Aipk to be 
symmetric, we finally arrive at 



Ph{n) 



1 1 



;cos2A^) n , 



with 



cos2A(p/ £ ) = / p(Ay)fe) cos 2Atpi s ,dAipi t 



(45) 



(46) 



If, for simplicity, we assume that the rotation angle can 
take just two values ±Aip, each with the probability 1/2, 
then Eq. ( [fBl) reduces to 



Ph{n) = i + icos"2A(^. 



(47) 



Equations (|45|) and ( |47| ) describe a purely exponential 
Markovian decay, for which we cannot expect any inhi- 
bition by the QZE. 



B. Simple non-Markovian model 

1. Description of the model 

Let us now consider a simple model for discrete non- 
Markovian dephasing, in which the probability of the ro- 
tation angle in the nth step depends on the rotation angle 
in the (n — l )th step (the so-called "random walk with 
persistence" 28 1). 

We assume that the rotation step can take two values 
±Aip with equal probabilities po(±Aip) = 1/2, whereas 
the nth rotation angle is equal or opposite to the pre- 
vious one (n — 1) with the probability p or q = 1 — p, 
respectively. Correspondingly, the conditional probabili- 
ties P(Aip n+ i, Aip n ) are denned by 



P(A<p n , A(p n ) = p, P(-Aip n , A(p n ) = q. (48) 

As a consequence, for a given A(p n the conditional aver- 
age of Aip n+1 is 



(A<f n+1 ) AiPn = ryAipn, 
j = 2p-l. 



(49a) 
(49b) 



The rotation steps are correlated for 1/2 < p < 1 (7 > 0), 
anticorrelated for < p < 1/2 (7 < 0), and uncorrelated 
for p = 1/2 (7 = 0). Equation ( t49a|) implies that the 
correlation function is given by Eq. (|28|) with 



B = Aip. 



(50) 



For small jumps and large n, Eq. (pOj), the evolution 
of the probability of horizontal polarization can be ap- 
proximated by combining Eqs. j2l|), (|3l|), and (|5C|). 
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2. Free evolution 



In the case of free evolution, Eqs. (C5) and (C6) in 
Appendix ^| yield the analytical solution 

u I \ 1 _l ff( r ) - 9(~r) 

where g(r) = (qcos2A(p + r)(p cos2Aip + r) n and r = 
\/ q 2 — p 2 sin 2 2Aip. 

We now consider special cases of Eq. (|5l|). If p = 
q = 1/2 we get the Markovian exponential decay given 
by Eq. p7|). For 4Atp 2 <C l,q 2 /p 2 one obtains r sa 
q — 2p 2 Aip r /q « g and Eq. ( |5l| ) reduces t o Eq ( p2| ) with 
i?o = .p(A</3) 2 /<7T r , which, in view of Eqs. ( 49b ) and (50), 
reduces to the formula for R in ( |3l| ) with (9=1. The 
corresponding behavior can be interpreted as a random 
sequence of "independent" rotations of the length A(pp/q, 
each rotation taking place after p/q steps. Using the same 
arguments as for the derivation of Eq. (E7p, we arrive at 



Ph(n) 



/ P 
cos 2Aip- 

V 9 



(52) 



Finally, for p = 1 (i.e., 7=1) the rotation steps are 
infinitely long - the exponential approximation is never 
true and Eq. (51) yields the "Rabi oscillations" (0). 



3. "Impulsive" (projective) measurements. 

Let us now assume that during each round-trip the 
polarization state of the photon is projected on the hor- 
izontal polarization (9 — in (||)). The total probability 
of the horizontal polarization is given by Eq. (||), inde- 
pendently of the memory of the rotator. 

Since for p > q this exponential decay is slower than 
the non-Markovian decay ( ^2|) [or, more exactly, (gj])], 
there exists a region of n values for which the probability 
of the horizontal polarization is larger with the projection 
than without it. (The region is finite because ( |52] ) and 
( |5i"| ) decay towards 1/2 whereas (||) decays to zero.) 

Figures ^(a)-(c) represent the discussed model. Figure 
|(a) exhibits the QZE. In Fig. |(b) we see the case of ex- 
ponential Markovian decay (p = 1/2) for which the QZE 
does not take place. It differs from the uninterrupted 
evolution, which tends to equilibrium of the population. 
Figure |(c) exhibits the AZE for p < 1/2. 

Generally, the smaller Aip the slower are thermaliza- 
tion and dec ay. For 1/2 < P < 1 (corresponding to 



< 7 < 1 in ( ]49b| ); cf. Sec. |V A| ) the larger p the larger 
is decay inhibition by the measurements (QZE), whereas 
for < p < 1/2 (—1 < 7 < 0) the opposite trend, i.e., 
decay acceleration (AZE) occurs. The generic properties 
of the QZE and the AZE clearly follow from this very 
simple model. 
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FIG. 4. Evolution of the probability of horizontal polariza- 
tion for A<p = 4°: (a) p = 0.8, (b) p = 0.5, (c) p = 0.3, as 
a function of the number of round trips n. Inset (a): same 
parameters as (a) but for smaller n. Inset (c): same param- 
eters as (c) but for smaller n. Solid line, the probability of 
the horizontal polarization without measurements, Eq. (pl[); 
dash-dotted line, the exponential approximation ([H^); dashed 
line, decay with perfect measurements (complete absorption), 
as in Eq. (||). 



VII. CONCLUSIONS 

We have obtained the general conditions of decay in- 
hibition or acceleration via polarization measurements. 
This has been accomplished by expressing the decay rate 



8 



as a convolution of two functions: (i) the fluctuation 
spectrum of the random polarization-rotation rate, which 
is analogous to the spectrum of the infinite-temperature 
reservoir; (ii) the Fourier transform of the measurement- 
induced polarization dephasing, which is determined by 
the absorber transparency. Inhibition of the decay by 
frequent measurements (the QZE) has been shown to oc- 
cur when the reservoir spectrum is peaked around zero 
frequency, which corresponds to correlated polarization- 
angle jumps. However, if the jumps are anti-correlated, 
the reservoir spectrum is split into two peaks, with a dip 
at zero frequency. Then one should obtain the quantum 
AZE, i.e., an increase of the decay rate as the rate of 
effective measurements increases. 

In [fl2"t only direct measurement schemes were consid- 
ered (Cook's U scheme with pulsed or cw measuring 
field), whereas here we have considered an indirect de- 
structive measurement scheme of the polarization state. 
Our analysis confirms that the general formula in fl2] | 
holds for direct and indirect measurement schemes, irre- 
spective of their differences, provided they approximately 
yield projections on the measured state. The general for- 
mula in |l2j ] for the effect of frequent measurements on 
decay was claimed to hold not just for ideal measure- 
ments, i.e., instantaneous projections, but also for con- 
tinuous measurements. Generally, measurements can be 
ineffective (i.e., producing projections with a probabil- 
ity less than one) ^3|. Here we have considered a com- 
prehensive measurement model, which encompasses the 
various types of measurements: ideal, ineffective (impul- 
sive), and continuous. The present analysis corroborates 
the generality of the formula in fl2|| and provides unify- 
ing expressions for the measurement-induced broadening 
function F(u>) [Eq. (^5|)j and the effective measurement 
rate v [Eq. @]. 

The significance of the present analysis lies in the phys- 
ical simplicity of the model and its experimental realiz- 
ability, as well as in the ability to control the reservoir 
spectrum (or memory time) and measurement rate by ad- 
justment of the Pockels-cell modulation (on a ns scale) 
and the PBS transparency, respectively. It allows us to 
test theoretically (and, hopefully, experimentally in the 
future) the general conclusions from Ref. [|l2| for the Zeno 
and anti-Zcno effects. 
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APPENDIX A: DERIVATION OF THE MASTER 
EQUATIONS AND THEIR VALIDITY 
CONDITIONS 

To derive the equation for the average solution of Eq. 
(0) we use the cumulant expansion technique p5],p6| • As- 



suming, for simplicity, that the odd moments of fi(t) van- 
ish, we obtain 



dQ 
dt 



= A(t)Q 



dt'K(t,t')Q(t'), 



(Al) 



where 



K(t,t') = e 2 {t,t')+ dh dt 2 e 4 {t 7 t 1 ,t2,t') + . 



(A2) 



Here the totally ordered cumulants are 



e 2 (i, t') = k(t, t')CU(t, t')C, (A3a) 
6 4 (i, h, t 2 ,t') = M 4 (t, ti, ta, t')CU(t, ti)CU(ti,t 2 ) 
xCU(t 2 ,t')C (A3b) 

where U(t, t') — exp^, dTA(r)] and 

M A {t,t x ,h,t') = {n{t)n{ti)^{t 2 )^{t')) - k(t,ti)k(t 2 ,t') 

(A4) 

is the the fourth cumulant of fJ,(t), the angular brackets 
denoting the average. 

Truncating the cumulant expansion ( |A2j ) at the sec- 
ond order, one obtains the following integro-differential 
equations for the elements of the polarization tensor, 

% = -2 f dt'k(t,t')0(t,t')[P h (t') - P v (t% (A5a) 



dt 



dl 



2 / dt'k(t,t')9(t,t')[P h (t') - P v (t')} 
Jo 

-T(t)P v , (A5b) 

-^-u-2 f dt l k{t,t')[l + 6 2 {t 1 t')]u{t l ). (A5c) 
2 Jo 



Note that Eqs. (A5a) and ( A5b| ) for the polarization 
probabilities are uncoupled from Eq. (A5c) for the co- 
herence. 

We assume that Ph(t') and P v (t') are slowly changing 
with respect to the integral kerne l. Th en on e can set 
Phv{t') « P hv (t) in Eqs. jA5a| ) and ( |A5b| ), yielding 
Eqs. (M). 



Consider now the validity conditions of Eqs. (10) for 
the case when fi(t) is a continuous-time stationary ran- 
dom process and T(i) = Tq is constant. First, the condi- 
tion to transform Eqs. (A5a) and (A5b) into Eqs. ( |To| ) 
can be seen to be 



r . 



(A6) 



Second, one can show that taking into account the fourth 
cumulant (A3b) in Eqs. (Al) and (A2) amounts to the 
addition of the term 

-A&dtx J* 1 dhMiitMMit'W^O^t 1 ) 

^[l + e^hM)] (A7) 







to the kernel 2k(t,t')0(t,t') of Eqs. jA5aD and ( |A5b| ). 
The cumulant M4(t,ti,t2,f) tends to zero with the in- 
crease of the difference between any two of its arguments. 
Assuming that all the characteristic decay rates are of the 
order of Tr, one can estimate that the fourth cumulant 
contributes to the rate R in Eqs. ([l8]) a quantity of the 
order of i? 2 /rR. Hence, the fourth cumulant can be ne- 
glected under condition (14). Note that condition ([u]) is 
stricter than flAq). 



n — 1 rn 



W„ = nB 2 + 2^ K " 



m=l m'=l 



In the case (28), the calculation of ( 
under condition ( |30a| ), 

I+76I 2B 2 -f6 



W n sa nB- 



1-70 (1-70) 5 



For 



(B5) 



shows that, 



(B6) 



APPENDIX B: GENERAL ANALYSIS OF SMALL 
DISCRETE PHASE JUMPS 

The correlation time T^ 1 of a random chain Aip n can 
be defined generally as 



Assuming the validity of condition ( 
Eqs. ( [To|) with the initial conditions ([ 
at t = tit-c in the form 



(Bl) 

|4f) , the solution of 
iOI) can be obtained 



Ph{t) 



= exp 



-W n - T t 



(B2) 



where W n = fitfRW). 

To prove Eq. (B2) we cast Eqs. 
q = L(t)q, where 



( [Toj ) in the form 



lit) 



Ph(t) 

Pvit) 



L(t) 



( -R(t) R(t) 

\ R(t) -R(t) - r(t) 



(B3) 



Hence q(t) = U(t)q(0), where U(t) = Texp[ J* dt'L(t')], 
T being the time-ordering operator. Note that L(t) ps 
L(oo) = L for t 3> T^ 1 . Consider two cases. Let, first, 
r <C Tr . We want to show that T can be omitted in the 
expression for U(t), resul ting in the coincidence of U(t) 
with the matrix in Eq. |B2j). For t < (R + r ) _1 one 
gets U(t) « [1 + j Q T dt'L(t')} ps eKp[j*dt'L(t% whereas 
for t > (i?+ro) _1 , choosing time r such that T^ 1 «C t <C 
(i^ + To)- 1 , one obtains U{t) Pa [1 + /J" dt'i(t / )]e i(t " r) « 
exp[J* di'L(f')]. The second case r > Tr implies, in 
view of (Jul), that To 3> -R. Then one can set P v (t) pa 
in Eqs. (pf), yielding i\ (t) pa e~ Wn , which can be shown 



to follow also from (B2). This finishes the proof of (JB2|). 

The above expression for W n , with the account of Eqs. 
(|j"l|)-(|l3|) and (||), can be cast as 



K„ 



\ m — rn ! \ 



(B4) 



or, for a stationary process, 



B 2 « (1 - 7 0) 2 



(B7) 



the second term in the right-hand side of Eq. (B6) can be 
neglected and Eq. (B2) yields (||l]) with the definitions 

For an arbitrary co rrela tion function K n one can show 
with the help of Eq. @ that Eq. (0) for R should be 
substituted by 



R 



(B8) 



or, equivalently, by Eq. (|32]). 

The validity conditions of the above res ults are: t ^> 
F^ 1 [roughly corresponding to inequality (30a)], the in- 
equalities (Jl4| ) and B 2 -C 1. In the case fl28|) the latte r 
two inequalities can be combin ed t o obtain Eq. ( 301 ), 
which is generally stricter than flBTj). 

Equation ( |32"|) can be rewritten in the form (|33[ ) with 



FW = I £ *H 



2tt 



(B9) 



n— — 00 



and 



G(w) 



1 



27TT r 
1 

2WT T 



n— — 00 



B 2 



cos riwr r 



(BIO) 



The function F(u)) (B9) is normalized to one in the in- 
terval — 7r/r r < uj < 7r / r r . Performing the summation in 
Eq. |g) yields Eq. @. In the special case @ Eq. 
( B1C ) reduces to the closed formula (|3) . 



APPENDIX C: EXACT SOLUTION FOR 
DISCRETE NON-MARKOVIAN PHASE JUMPS 

Here we calculate polarization dephasing for a rather 
general model of random rotation angles. The average 
phase factor in Eq. ( |43| ) is given generally by 

( e «vn) = e 2 ^"p(A^,...,A^ n ), (CI) 
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where p(Aipi, . . . , Atp n ) is the joint probability of the 

random chain Atpk- 

According to the model of Sec. VI B 1 , tp n — J2k=i ^^fc 
is a (non-Mar koyian) sum over a Markovian chain Aip^- 
The average (Gl) can be calculated for arbitrary Marko- 
vian chains, i.e., chains in which Aifk+i depends only on 
Acpk, as follows. In this case 



p(A<fii, . . . , Aip n ) = P(A(p n , Alfin^x) . . . 

x?(A^,A^)p (A W ; 



(C2) 



Here po(Aipi) is the unconditional probability and 
P(A(pk, Aifk—i) is a conditional probability. Combining 
Eqs. flC^), (|Cl|) and (ph, one can obtain that 



(C3) 



The quantity / n +i(A(/5„ + i) is defined by the following 
iterative relation, which is convenient for numerical cal- 
culations, 

f 1 (Acp 1 )=p (A<p 1 ), 

f n+1 (A^ n+1 ) = ]T P(A<p n+1 , A^)e 2lA ^/„(A^). (C4) 

The above solution can be conveniently written in a 
matrix form. Assuming that the random variable Aipk 
can assume the values 6(pi (i — 1,2,...), we represent 
Po(A(pk) by a column vector pg with the components 
Po(Sifi). The conditional probability is represente d by 



the matrix P with P i3 = P(8tpi,5<pj). Then Eq. (C3) 



and (|CJ) yield 

Ph(n) 



-Re[u(P$)> ] 



(C5) 



where u is the row vector with the components equal to 
1 and $ is a diagonal matrix with (f>j, = e 2l5 ^ : > . 



For the non-Markovian model of Sec. VI B 1 



Po 



P = 



P 1 
q p 



$ = 



e 





-2iAip 
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